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ABSTRACT 

We evaluate the large scale polarization of the cosmic microwave background induced, via 
Thomson scattering prior to decoupling, by nearly scale- invariant spectra of scalar and tensor 
metric perturbations, such as those predicted by most inflationary models of the early Uni- 
verse. We solve the radiative transfer equation for the polarized photon distribution function 
analytically, for wavelengths of order or larger than the horizon size at decoupling, assuming 
no reionization. The induced polarization is proportional to the redshift induced by the per- 
turbations at decoupling and to the duration of the decoupling transition. Normalizing the 
spectra to the quadrupole anisotropy measured by the COBE satellite, the expected degree 
of linear polarization is P < 10^^, more than two orders of magnitude smaller than current 
bounds. The dependence of P with the tilt in the spectrum away from scale invariance is 
too weak to be of observational relevance or to provide a sensitive discrimination between 
scalar and tensor fluctuations. The variation of P with the angular scale is significant. 
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The most plausible explanation of the quadrupole anisotropy in the temperature of the 
cosmic microwave background measured by the COBE satellite is that it originates in 
redshifts induced by metric fluctuations around a Robertson- Walker background, through 
the Sachs- Wolfe effect ||. 

The fluctuations in the metric may be scalar (energy-density) or tensor perturbations 
(gravitational waves). Inflationary cosmological models predict them both, as a result of 
quantum fluctuations at early stages. Exponential inflation predicts a scale-invariant, gaus- 
sian spectrum of scalar fluctuations 0, and a smaller amount of tensor fluctuations [Q. 
Other inflationary models, for instance power-law inflation , predict spectra slightly tilted 
away from scale invariance. It would be of great value, and would shed much light upon the 
Universe early evolution, to know what fraction of the quadrupole anisotropy measured by 
COBE is due to scalar fluctuations and what fraction is due to gravitational waves. Most 
inflationary models predict a ratio between scalar and tensor fluctuations proportional to the 
tilt in the spectrum away from scale invariance. In principle, independent measurements of 
the relative amplitudes of scalar and tensor fluctuations and of the spectral index, could thus 
provide significant evidence for inflation |^, H, |^, if they happen to fit the predicted 
ratio. Put it another way, independent measurements of these parameters would allow a 



determination of the potential of the field driving inflation This program however, is 
not easy to carry out in practice, among other things because of the statistical nature of 
the predictions (the "cosmic variance" problem), and because discrimination between tensor 
and scalar sources of anisotropy requires measurements at small angular scales |]^, [1^ , where 
the primordial spectrum is affected by many different processes that strongly depend on the 
details of the cosmological model. 

Its degree of polarization is another (as yet unmeasured) property of the microwave 
background radiation, which determination would shed more light upon the Universe early 
evolution. The polarization of the microwave background radiation is actually tightly re- 
lated to its anisotropy: Thomson scattering of anisotropic radiation by free electrons before 



recombination induces some degree of linear polarization [|^, Q 



In this paper we evaluate the degree of linear polarization induced both by scalar as 
well as tensor fluctuations, with a gaussian, power law spectrum, not far away from the 
scale invariant (Harrison-Zeldovich) type. Our analysis is based on the analytic method 
and approximations developed by Basko and Polnarev to study the polarization induced 
by anisotropic expansion of the Universe [|T3| and by gravitational waves |]TB|]. We extend 



their method to include the effect of scalar fluctuations. Closely following their approach to 
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the subject, we solve the radiation transfer equations for the polarized photon distribution 
function analytically, performing approximations which will prove valid for the large angular 
scale (> 1°) polarization. We assume that the Universe is already matter-dominated at 
recombination, and that there is no further reionization. Polarization measurements of the 
cosmic microwave background at large angular scales could very well serve as a test of the 
reionization history [|14|, |16|, [1^, but this is not the focus of the present paper. Instead, we 



concentrate (assuming no reionization) on the predictions of inflationary models, such as 
power law inflation. In other words, we normalize the scalar and tensor fluctuations spectra 
to the COBE measurements of the anisotropy in the microwave background temperature, 
and determine the expected large angular scale polarization. From Polnarev's results ||15[ 
one expects a rough estimate on the degree of polarization of the microwave background 
of around a few percent of its large scale anisotropy, and this is a rather small number, 
about two orders of magnitude smaller than current bounds. We feel, however, that the 
recent measurement of the quadrupole anisotropy justifies a more detailed analysis of the 
inflationary predictions for the large angular scale polarization. Our motivation is to find the 
dependence of the large scale polarization with the spectral index of the fluctuations, which 
if large could in principle provide an adittional indirect test for inflation. Other dependences 
of the degree of polarization, such as upon the angular scale or the cosmological parameters, 
are also of interest. 

Our work partially overlaps with a recent preprint by Crittenden, Davis and Steinhardt 
p!7| , where the same subject was treated numerically. Our approach, valid only for large 
angular scales under the assumption of no reionization, is of a more limited applicability 
than their full numerical treatment. On the other hand, our analytic results, based on the 
application of Basko and Polnarev's method [|1^, |15|, provide some interesting insights into 



the problem, such as the explicit dependence upon the spectral index, the angular scale and 
upon parameters of the cosmological model. 

The polarization of the cosmic microwave background radiation will be determined from 
the Stokes parameters that result from solving the Boltzmann equation in a perturbed 
spatially-flat Robertson- Walker background. The metric reads 

ds"^ = a'^{T)[-dT^ + {6ij + hij)dx'dx^] , (1) 

where r is the conformal time, related to proper time t hj dr = dt/a{t), \hij\ << 1, and 
i,j = 1,2,3. The polarized photon distribution function is represented, in a general case, 
by a column vector n(z/, 9, (p) = {Ii, Ir, U, V), its components being the Stokes parameters of 
the radiation that arrives from a direction defined by the angles 6,(f) with frequency u [|18|. 
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Ii , Ir are the intensities in two orthogonal directions defined by the unit vectors 6 and of 
spherical coordinates. The Boltzmann equation reads 

I rl j-2-K 



d i d \^ du dh ^ 

dr ^ dx"^ I dr du ^ ^ 



(2) 



Here P{fi,(j), fi',(j)') is a 4 x 4 matrix that characterizes the effect of Thomson scattering 
on polarization [0, fi = arccos^, e* is a unit vector in the direction of {0,(f)), ax is the 
Thomson scattering cross section and Ne is the free electron density. The effect of the metric 
perturbations on the frequency of the photons is given by the Sachs- Wolfe formula, which in 
the synchronous gauge reads [0 

Idu _ IdK^ 
ud^~2^''^ 

The Stokes parameter V actually decouples from the rest, and we will not take it into further 
consideration, since it is not necessary to determine the degree of linear polarization, given 
by P = VQ' + Uyi, with I = Ii + I^^Q = Ii- J,. 

We describe the growing mode of adiabatic scalar fluctuations in a matter-dominated, 
spatially-flat Robertson- Walker background, in terms of Bardeen's gauge-invariant variable 

( coincides with the energy-density fluctuations Sp/p in a matter-dominated regime for 
wavelengths inside the horizon. We decompose it in plane waves, ({x) = J d^ke'^'''^({k), where 
({k) is a gaussian, random variable with statistical average {C{k)({q)) = Pi^{k)6 {k — q) / Airk^ . 
deflnes the spectrum of scalar fluctuations, which we will take to be of the form 

P^k) = P^k-~' (5) 

where P^ is a constant, and the spectrum is scale invariant (constant amplitude on scales 
equal to the horizon at any given time) if n = 1. 

Analogously, tensor fluctuations in the transverse traceless gauge, that enter the horizon 
during a matter-dominated regime, are described by 

h^A^, ^)=T.J dV-'h,{k)e,^{X, k) (^^^^ (6) 

where A = 1, 2 denotes the gravitational wave polarization, characterized by e^j, ji are spheri- 
cal Bessel functions, and hx are random variables with statistical average {hx{k, T)hx' {q, r)) = 
Ph{k)S{k — q)6x\'/4:TTk^. We will consider a power spectrum of the form 

Ph{k) = Phk^-' , (7) 
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which is scale invariant if n = 1. 



We first discuss the polarization induced by one single scalar mode, of the form = 
i^{k)e'^^'^ . Since the problem has rotational symmetry around the direction of /c, only two 
Stokes parameters need to be considered {U = 0). Now fi is a two-component column 
vector, which we expand to first order in the perturbation, as fi = no + (l/2)noni(a;, r, z/), 
with fio = (l/2)no(l 1), where no(i^) is the blackbody distribution. Boltzmann's equation 
(0) simplifies considerably: 

2(l-/i2)(l_^'2)+^2^'2 ^2 



where q = axN^a, p = ^j^^ ^ 1 in the Rayleigh- Jeans zone of the blackbody spectrum, 
and Fs = —T(^{k)e^^'^k'^ /lb arises from the redshift induced by the perturbation, as given 
by the Sachs- Wolfe equation (^). We will solve equation (|^) in the long wavelength limit, 
neglecting the spatial derivative in the l.h.s., but anyhow keeping the k ■ x-dependence in 
the distribution function. This amounts to a quasi-stationary approximation []T3|, valid for 
perturbations with wavelength much longer than the photon mean-free path just prior to 
decoupling {k « q), and also much longer than the distance that a photon can travel 
during the decoupling transition, as we will discuss in more detail later. We normalize the 
present conformal time Tq = 1, then a(r) = and q ~ 0.14i7;,/iXe(l + zY, with the 
density in baryons in units of the critical energy density, Xf, the ionization fraction, and 
Hq = 100 h km seg~^ Mpc"^. Thus around decoupling, q ^ 10^. Notice that with our 
normalization k 271 corresponds to a wavelength comparable to the present horizon, and 
k ^ 30(27r) to the horizon at decoupling. 

We solve eq. @ through the following ansatz for ni: 

fii = a(r)a + /3(r)b + 7(r)c (9) 

where 

a=(/i^-^)Q ; b = (i-^2)(^_^J ; c=Q (10) 

With this ansatz, and neglecting the spatial derivatives in the l.h.s. of the Boltzmann 
equation (H), we derive the following equations 

9^ + ^^ = ^^ ' 9^+10^^ = -10^ ' d-r=-/' ^''^ 

where we have defined the variable ^ = a + (3, which decouples from j3, and is a measure of 
the anisotropy in the photon distribution function. f3 instead is a measure of the polarization. 
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since Ii — Ir = noP{l — /i^), and the anisotropy, as measured by ^, is its source. 7, on the 
other hand, decouples from the rest, and just reflects a local adjustment of the black body 
intensity, independent of the direction of propagation. The equations for ^ and P can be 



formally integrated as follows [15 



e(r) = / Fs(r')e-''(^'-')rfr' (12) 



^(^) = -T^r q{r')aT')e-^Ary)dr' (13) 
lU JO 

where k{t,t') = J^, q{T")dT" is the optical depth. We wish to evaluate /3(ro) from eq. ([T3|) . 
This equation shows that the present polarization of the microwave background (assuming 
no reionization) is essentially that produced at the times when qe~^'^^'^°''^''> is significantly 
different from zero, i.e. around the time of recombination, since much later the free electron 
density (and thus q) is negligible, while much earlier the optical depth is very large. To be 
more precise, the present polarization is the result of Thomson scattering around the time of 
decoupling of matter and radiation, which occurs slightly after the free electron density starts 
to drop significantly. Although these values depend on the particular cosmological model 
considered (through Q, Qb, Hq, etc.), typically recombination occurs around z ~ 1400 while 
decoupling takes place between z ~ 1200 and z ^ 900. This is the interval during which the 
differential visibility differs significantly from zero, and is actually well approximated by a 
gaussian [^]. We will denote by the conformal time around which decoupling occurs, and 
Ard the duration of the process. These quantities actually do not depend very strongly on 
the cosmological parameters, typical values being (with the convention tq = 1), ~ 3 x 10~^ 



and ATd^Sx IQ-^ |20, 21 



As explained above, we only need to know ^(r) around in order to integrate eq. ([T3|) . 
If we define k{t) = k{tq,t), then k{t,t') = k{t') — k{t), and we may approximate around 

e(r) ^ Fs{T,)e<-^ r e-^^^'Mr' ^ Fs{Ta)AT,e^^^~^ E{k{t)) (14) 
Jo 

where in the last step we made the approximation "^^ich is justified by the 

almost gaussian shape of the differential visibility around decoupling, and we defined E{h) = 
Ji° ^~^dx. Then, using that —qdr = dK, we integrate /3: 

1 3 
P ^ -—Fs{Td)ArJ e-To-e^E{K)dK^-0.l7FsiTd)ATd . (15) 
10 JO 

Finally, the polarization produced by a single scalar mode is given by 

pi = mi-f^r ■ (16) 



6 



The polarization induced by a single tensor mode of a given polarization can be calculated 
in a similar way. Now, as opposed to the scalar case, there is no longer rotational invariance 
around the direction of k, thus U ^ 0, and n is a three column vector {Ii, Ir, U), which we 
expand as n = no + (l/2)nofii(x, r, z/, 0), with fio = (l/2)rio(l 1 0). The angular dependence 
of the redshift produced by the gravitational wave due to the Sachs- Wolfe effect suggests an 
ansatz for ni of the form: |jl5 

ni = a(r)a + /3(r)b (17) 



where 



a = (1 - fi"^) cos(20) 



/1\ 
1 



/ (1 + fi^) cosi2(j)) \ 

-(l + /i2)cos(20) 
V 4yusin(20) / 



This corresponds to a gravitational wave with polarization h+. The formulae for a polariza- 
tion /ix are identical to these after interchange of cos(20) with sin(20). With this ansatz, 
the equations for ( = a + (3 and /3 are exactly the same as those in eq.(|Tl|) for the scalar 
case, but with Fs replaced by Fnp = ^e^'''^hx{k) Under the same approximations 

assumed valid for the scalar case, the result for /? is then /? ^ 0.17Fr(rrf)Ar^. Finally, the 
polarization induced by a tensor mode is given by 

pi = i/3n(i+/.Y+V] (19) 

where we have already summed over the two polarizations of the gravitational wave. 

Both for scalar as well as tensor fluctuations, /3 oc F^TdjArd, with the appropriate F 
in each case. This result for (3, which is a measure of the induced polarization, has an 
intuitive explanation, after noticing that F measures the anisotropy in the rate of change in 
the frequency of the cosmic background photons. Indeed, the Sachs- Wolfe formula for the 
scalar and tensor modes can be written as 

-— = Fs^i^ ; -— = Ft 1-/.' COS 20 . 20 
V dr V dr 

This anisotropy in the redshift translates into a difference in the number of photons of 
a given frequency traveling in different directions, and this is what makes possible a net 
polarization after Thomson scattering. Besides, the polarization builds up over a period 
Ar^ only, since anisotropics in the photon distribution produced much before decoupling 
are erased by successive scatterings, while after decoupling scattering is negligible and thus 
polarization can no longer be induced. 

Next we evaluate the rms total polarization that results from the superposition of all the 
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statistically independent contributions of scalar and tensor modes respectively: 

{P^)^ = fd'k[{l + f^'y + A^']\P\' = lJ,'-^f{P,k--')[0A7Ar,£C^ 

We have introduced a cut-off, kmax, determined by the minimum wavelength to which a 
given experimental setup is sensitive. In general, measurements result from an average over 
a certain angular scale, for instance the aperture of the antenna's horn. With our conventions, 
a wavelength with comoving wavevector k subtends an angle 6 such that k ^ tt/ sin(0/2). 
This relation determines kmax for measurements averaged over an angle 6. 

Finally, for scalar perturbations 

^ ^ 15 (n + 3)^ 15 ^ ^ ^ 

while the result for tensor perturbations can be written as 



{P^)t = Dpin,kmax){P^)s , (23) 



where 
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Dp = 13.5 g{n,kmax)^ ; ^mao;) = 25(n + 3) — ^3 ■ (24) 

With this definition, g{n, kmax) 1 when kmax^d ^ 1, ^-e. for wavelengths well outside the 
horizon at decoupling. The total rms polarization is given by the sum of the contributions 
due to scalar and tensor modes, given their statistical independence: 



^) = .^{P^)s+{P^)t (25) 



Now we assume that the scalar and tensor perturbations have as common origin a period 
of power law inflation, in which case the scalar and tensor power spectra are related by 
P/i/P^ ~ 4(n — l)/9(?2 — 3) 0, H, 1^. The ratio of the tensor and scalar contribution to the 
quadrupole moment in the temperature anisotropy can also be approximated by a simple 
function of n, {af) t / {o-f) s = Da ^ 13{n — l)/(n — 3) 0. Assuming that the anisotropy 
measured by COBE is due to these inflation-produced perturbations alone, then it is simple 
to express both P,^ and Ph in terms of the measured quadrupole and the spectral index. For 
instance: 

^ ^ 2^-"(3-n) r((9-n)/2) r^(2-n/2) 
^ 47r2 (8-7n) r((3 + n)/2) r(3 - n) ^ l^^j 
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with {a\) = {al)s + {a\)T ~ 2 x 10~^, the quadrupole measured by COBE. Finally, the total 
rms polarization predicted by power-law inflationary models, as a function of the spectral 
index n, is given by 

P - 1.2 X 10-^ ( \ ( ^1 ' ( \ i Cin, A:_) (27) 

V3 X 10-3/ V 50 ; V3 X 10-2/ ^2 x 10-^ J ^ ' """"^^ ^ ' 

with 

« ^„„ir((9-n)/2)r2(2-n/2) fl + D, 



3(r2 + 3)7r '"""2"-ir((3 + n)/2) r(3-n) U + Z^a/ ^ ' 

The function C has been defined such that C = 1 if n = 1. The value kmax ~ 50 corresponds 
to measurements that average over an angular scale of order Q ^ 7°. 

In Fig. 1 we plot the total rms polarization for k^ax = 50, as well as the individual tensor 
and scalar contributions, for values of n between 0.5 and 1, those consistent with COBE 
measurements. The predicted degree of polarization is more than two orders of magnitude 
smaller than current bounds, P < 6 x 10~^ [22|, obtained with a 7° aperture horn at 33 GHz. 



The polarization induced by tensor modes vanishes if n = 1 and equals the scalar contribution 
when n = 0.5. The total polarization decreases by a factor about 0.6 as n varies from 1 to 
0.8 {n ^ 0.8 corresponds to equal contributions of scalar and tensor modes to the large scale 
anisotropy), and about 0.4 between n = 1 and n = 0.5. Notice that the dependence of 
the total polarization upon n is mainly due to the tilt in the spectra, i.e. it comes mostly 
from the term kf;^~J-^^'^. This is because polarization is proportional to the redshift induced 
by the perturbations at decoupling, and is no longer produced afterwards, while anisotropy 
continues to build up. Consequently, the dominant wavelengths for polarization are much 
smaller than those that dominate the large scale anisotropy. Once the spectra are normalized 
to the large scale anisotropy measured by COBE, smaller n implies less power on smaller 
scales. 

A large class of inflationary models predict a definite ratio between the spectral index 
n and the scalar and tensor contributions to the anisotropy in the microwave background 
1^, . We have shown, applying the methods developed by Basko and Polnarev [1^ to 



the specific spectra predicted by inflationary models, that there is also a definite relation 
between the net polarization and n, at least assuming no reionization. Unfortunately, the 
rather weak dependence of P with n makes it unlikely that measurements of polarization 
could provide a sensitive confirmation of these predictions. Our conclusions in this regard 
are similar to those of Crittenden, Davis and Steinhardt [0, that in a recent preprint 



approached the same subject numerically. Our analytic results, on the other hand, are only 
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valid if kmax^Td << I7 which means kmax << 300, while the numerical treatment of Ref. 
IJT^ applies to a more general situation. 

The large scale polarization depends upon cosmological parameters, such as fl, Qb, Ho, 
etc., only through Td and Ar^, i.e. through the parameters that characterize the decoupling 
transition, and this is a very weak dependence . 

Notice finally that P is proportional to k'^ax- This is a consequence of the fact that the 
polarization induced by a single mode is proportional to the redshift it causes at decoupling, 
irrespective of its wavelength, and because in the linear regime the perturbations grow as 
(fcr)^. Thus, the dependence of the degree of polarization with the angular scale over which 
polarization data is averaged is significant. 

References 

[1] G.F. Smoot et al., Ap. J. 396 (1992) LI. 
[2] R. Sachs and A. Wolfe, Ap. J. 147 (1967) 73 

[3] A.H. Guth and S.Y. Pi, Phys. Rev. Lett. 49 (1982) 110; S.W. Hawking, Phys. Lett. 
B 115 (1982) 295; A.A. Starobinsky, Phys. Lett. B 117 (1982) 175; J.M. Bardeen, RJ. 
Steinhardt and M.S. Turner, Phys. Rev. D 28 (1983) 679. 

[4] V. A. Rubakov, M. Sazhin and A. Veryaskin, Phys. Lett. B 115 (1982) 189; R. Fabbri 
and M. Pollock, Phys. Lett. B 125 (1983) 445; L. Abbott and M. Wise, Nucl. Phys. 
B 244 (1984) 541; L. Abbott and D. Harari, Nucl. Phys. B 264 (1986) 487; B. Allen, 
Phys. Rev. D 37 (1988) 2078. 

[5] L. Abbott and M. Wise, Nucl. Phys. B 244 (1984) 541; F. Lucchin and S. Matarrese, 
Phys. Rev. D 32 (1985) 1316; R. Fabbri, F. Luchin and S. Matarrese, Phys. Lett. B 166 
(1986) 49. 

[6] L. Krauss and M. White, Phys. Rev. Lett. 69 (1992) 869. 

[7] R.L. Davis et al., Phys. Rev. Lett. 69 (1992) 1856. 

[8] F. Lucchin, S. Matarrese and S. Mollerach, Ap. J. 401 (1992) L49. 

[9] D. Salopek, Phys. Rev. Lett. 69 (1992) 3602. 



10 



[10] A. Liddle and D. Lyth, Phys. Lett B 291 (1992) 391. 

[11] E.J. Copeland, E.W. Kolb, A.R. Liddle and J.E. Lidsey, "Observing the infiaton poten- 
tial", preprint FNAL-Pub-93/071-A, (1993); M.S. Turner, "Recovering the inflationary 
potential", preprint FNAL-Pub-93/182-A (1993). 

[12] R. Crittenden, J.R. Bond, R.L. Davis, G. Efstathiou and P.J. Steinhardt, Phys. Rev. 
Lett. 71 (1993) 324. 

[13] M. Rees, Ap. J. 153 (1968) LI; G.P. Nanos, Ap. J. 232 (1979) 341; J. Negroponte and 
J. Silk, Phys. Rev. Lett. 44 (1980) 1433. 

[14] M. Basko and A. Polnarev, Sov. Astron. 24 (1979) 3. 

[15] A. Polnarev, Sov. Astron. 29 (1985) 607. 

[16] K.L. Ng and K.W. Ng, "Large scale polarization of the cosmic microwave background", 
preprint of the Institute of Physics, Academia Sinica, Taiwan, IP-ASTP-08-93 (1993). 

[17] R. Crittenden, R.L. Davis and P.J. Steinhardt, "Polarization of the microwave back- 
ground due to primordial gravitational waves". University of Pennsylvania preprint 
UPR-0575T ( |astro-ph/9306027D (1993). 

[18] S. Chandrasekar, Radiative Transfer, Oxford University Press (1950). 

[19] J. Bardeen, Phys. Rev. D 22 (1980) 1882. 

[20] J.R. Bond, in The Early Universe, edited by W.G. Unruh and G.W. Semenoff, D. Reidel 
Publishing Co. (1988). 

[21] R. Sunyaev and Y. Zeldovich, Astrophys. Space Sci., 7 (1970) 3. 

[22] P. M. Lubin and G.F. Smoot, Ap. J., 273 (1983) LI. 

Figure Caption 

Fig. 1: Total rms large scale (~ 7°) polarization of the cosmic microwave background 



P = y-Rf + -P| and individual scalar and tensor contributions, Ps and Pt, as a function of 
the spectral index n in power-law inflationary models, assuming no reionization. 
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